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Abstract 

Phase transitions generically occur in random matrix models as the param- 
eters in the joint probability distribution of the random variables are varied. 
They affect all main features of the theory and the interpretation of statistical 
models . In this paper a brief review of phase transitions in invariant ensembles 
is provided, with some comments to the singular values decomposition in com- 
plex non-hermitian ensembles. 



1 Phase transitions in invariant hermitian en- 
sembles. 

Random matrix ensembles have been extensively studied for several decades, 
since the early works of E.Wigner and F.Dyson, as effective mathematical ref- 
erence models for the descriptions of statistical properties of the spectra of 
complex physical systems. In the past twenty years new applications spurned 
a large literature both in theoretical physics and among mathematicians. Sev- 
eral monographs review different sides of the physics literature of the past few 
decades, such as@,@,@||||l|,||||@. Their combined bibli- 
ography, although very incomplete, exceeds a thousand papers. Sets of lecture 
notes are fl0§ §§ @@ |§ @ . The classic reference is Mehta's book |§. 

For a long time studies and applications of random matrix theory in large 
part were limited to the choice of gaussian random variables for the independent 
entries of the random matrix. This was due both to the dominant role of the 
normal distribution in probability theory as well as to the nice analytic results 
which were obtained. Increasingly, in the past two decades, a wide variety of 
matrix ensembles were considered, where the joint probability distribution for 
the random entries depends on a number of parameters. As the latter are al- 
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lowed to change, the generic occurrence of phase transitions emerged. 

In this section, I begin by recalling the problem in the easiest case, the invariant 

ensemble of hcrmitian matrices . 



Let H — (Hij)i j = t t „ t ff be hermitian random matrix with joint probability 
density for the independent entries 

P{Hu,H 12 ,..,H NN )dH = e- N ^nH) dH l J e -NTrv(H) dH 

V(H) = X -a 2 H 2 + \a A H± + .. + ^a 2p H 2 P , a 2p > , 
dH = YiidReHijdlmHij) Y[dH u (1.1) 

i<j i 

We are interested in the partition function Z^{a 2 ^ a A , .., a 2p ), the free energy 
Fn{o, 2 , 0,4, .., a 2 p) = — log Z^(a 2 , a4, .., a 2p ), the "one point resolvent" Gn{z) , 
the connected correlator G^{z\, z 2 ) 

Z w (a 2 ,a 4 ,..,a 2p ) = J e ~ NTlV ^ dH , G N (z) = Itt < -± 

Gff(zi,za) =< Tr-L-Tr-L- > - < Tr-L- >< Tr-L- > 

(1.2) 

The name invariant ensemble for the ensemble of these hermitian matrices re- 
minds that since the density P (Hn, Hi 2 , .., Hnn) is invariant under a similar- 
ity transformation H —> XJHU~ X with arbitrary unitary matrix U, most of the 
interesting quantities, like those in eqs.( |l.2| ) may be evaluated from the joint 
probability density of the eigenvalues. 

Also important are the monic polynomials P n {z) = z n + 0(z n ~ 1 ), orthogonal 
on the real line with the weight e~ NV ^ 



H > 



P n (z)P m (z)e- NV ^dz = h n S, 



n u nm i 



zP n (z) = P n+1 (z) + R n P n - 1 (z) , R n =p^ > Q , 

N—i 

Z N (a 2 ,a 4 ,..,a 2p ) = M (h ) N J] (fl*)" - " (1.3) 

n=l 

One obtains a non- linear recursion relation for the coefficients R n . For instance 
if V(x) = ^fx 2 + 2fx 4 + ^x 6 one has the recurrence relation (see |@ j H > 
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1 ) 



n 
N 



a 2 + a 4 (Rn-i + R n + R n +i) + a & (R n -i + R n + Rn+i) + 



(1.4) 



occasionally called " pre-string equation" or the Freud equation. One also intro- 
duces the set of orthonormal functions ij) n (z) and the two point kernel Kj\r(x, y) 
in terms of which all n— point correlation functions are expressible 



ipn(z) 



1 



N-l 



-NV 



(1.5) 



If all the coefficients a 2 k > the statistics of the eigenvalues may be evaluated 
in the limit N — ► oo , sec jl3J , limjv^oo Gn(z) = G(z) is holomorphic in the 
complex z plane, except for a segment (—A, A) on the real axis. Furthermore 



G(z) 



■A Z — \l 



p(jt) = Jim — Tr < S(p - H) >= Jim K N {p,p) 

N—>og JV N—>oa 



The limiting density of eigenvalues p(X) , the one point correlation function, is 
the unique solution of the integral equation : 



V'(X) = 2V 



d/z- 

-A A - /X 



p(A) dA 



(1.6) 



p(A) > on its support (—A, A) , it vanishes as a square root at the boundary 
p(X) ~ (A - IAD 1 / 2 . Further more the coefficients R n approach a smooth limit 
R{§) ~ R(x) and 



F(a 2 ,a 4 , ..,a 2p ) = iV 2 



dA p(A) y(A) 



dXdp P {X)p{p) log |A- M | + 



dx (1 - x) log i?(x) - 1 log fto + o(^; 



The free energy F(a,2, 04, .., a 2p ) is analytic in the couplings (a 2 , 04, .., a 2p ). The 
saddle point solution is equivalent to the resummation of the planar graphs, it 
is equivalent to the solution from the recurrence relations for R n and to other 
techniques, such as the loop equations. The orthogonal polynomial technique 
and the loop equations are superior to evaluate in a systematic way the terms 



in the series in the parameter 



(N 2 



corresponding to the resummation of the 



graphs which are embeddeble on orientable surfaces with k handles [101] 
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It was also proved that the connected two point correlators exhibit two 
different forms of universality (that is independence from the set of coefficients 
{ci2k} but depend only on the endpoints ±A) : 

Global universality. The limiting connected density-density correlator p c (A, A') , 
after smoothing over a scale much larger than the level spacing A^v , is [^0| 

n 1 A 2 -XX' 

Local universality. The limiting two point kernel K(x, y) has the sine law [ p0| 
for eigenvalues in the bulk of the spectrum, measured in units of A^r 

/ / \ sin[7r(s — s')] , . . . , . . A 

^buik(s, a!) = — V , s = X/A N , a! = X' A N 

tt(s — s ) 

or the Airy law , close to the tail of the spectrum (the soft edge) E3| || 

s s \ J\. 

A general derivation of spectral correlators is given in the recent paper by 
Kanzieper and Freilikher j70|. By following Shohat J97[ , the recurrence rela- 
tion for the orthogonal polynomials P n (x) is turned into an exact second-order 
differential equation for the orthogonal functions ip n (z) 

%^-^(A)^ + a„(AK(A) = 

Because of the smooth behavior of the coefficients R n , in the case where the 
eigenvalues have support on a single segment, the complicated forms T n {X) , 
Q n {X) simplify at large order and the differential equation leads to the global 
and the local universality results mentioned above. [] 

All the above is deeply affected by phase transitions. 

The limiting eigenvalue density p(X; 02, 04, .., a2 P ) continued to negative val- 
ues for one or several coefficients d2k (while keeping ai v > ) may not be 
positive definite. Then the integral equation ( |l.6| ) allows new solutions with 
eigenvalue density positive definite with support on two or more segments of 

If a logarithmic singularity is present at the origin, the Bessel law is derived 

, 7T ,,, J a + l/2(Ks)J a _ 1 / 2 (Trs') - J a - 1 /2(Ks)J a + 1 / 2 (7Ts') 

K oligin (s,s ) = - (ss yi* 

2 s — s' 

where s and s' are scaled by the level spacing Ajv(O) near the spectrum origin, s = A/Ajv(0) 
, and a > —1/2. For the simpler situations where the orthogonal polynomials are classical 
the Bessel law had been derived in H E7I H Hi . 
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the real axis. The correct solution minimizes the free energy. As one explores 
the space of the parameters, the free energy is evaluated on different saddle 
point solutions, which may coincide for certain critical values of the parameters. 
Then the free energy is usually a continuous but not analytic function of the 
parameters. 

The recursion coefficients R n no longer have a smooth "continuous" limit, which 
makes more difficult the orthogonal polynomials solution. 

The lack of analyticity of the free energy at critical values of the parameters 
is analogous to phase transitions related to spontaneous symmetry breaking in 
classical statistical mechanics. For instance, in the simplest case of potential 
V{x) = ^x 2 + CJ fx 4 , first analysed for negative values of a 2 in papers 96 |p5fl . 



it is convenient to add a linear term , which explicitly breaks the Zi symmetry, 
V{x) — a%x + ^x 2 + ^fx 4 , next perform the "thermodynamic limit" N — > oo, 
finally remove the symmetry breaking term a± — > 0, see J62[ , pfl . This allows 
the evaluation of the " order parameter" < Tr_ff > 

lim lim < TyH > = sign(ai) 0(—a,2 — lyfa^) f (122, 0S4) 

ai— >0jV— >oo 

For the simplest case V(x) = ^x 2 + ^fx 4 , it was shown that the correct ansatz 
for the recursion coefficients R n , if 02 < —2^fal is that the even i?2„ and 
odd i?2n+i approach two different smooth "continuous" functions pi. This 
period-two ansatz requires a little generalization in the case which includes the 
infinitesimal symmetry breaking term pij [^3| because it leads to recurrence 
relations 

zP n (z) = P n+1 (z) + S n P n (z)+R n P n - 1 (z) , R n = J^>0 , 
and it explains the origin of the period- two ansatz. 

However in the next simplest case , like V(x) — ^fx 2 + ^fx 4 + ^x 6 , or 
higher order polynomials, corresponding to multiple well potentials, which may 
not be degenerate, the behaviour of the coefficients R n is erratic and it is dif- 
ficult to reproduce the results of the saddle point analysis by the orthogonal 
polynomials || @@ || @. 

The lines of "phase transition" in the parameter space, related to the continu- 
ation to negative values of the coefficient d2 P of the monomial of highest order 
may be found in terms of previously discussed phase transitions by the addition 
of an infinitesimal monomial e x 2p+2 , see for instance [|7| |6^| . 

Connected correlators, when the support of the eigenvalue density is two 
segments, were shown || |l| to have a different form of global universality, 
involving elliptic integrals. In the case of multicritical behaviour the local uni- 
versality form has a modified Bessel law [|| . 
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Important recent works seem to be so powerful and comprehensive to solve 
the above mentioned ambiguities. The Freud equation is expressed in a matrix 
Lax representation and the semiclassical asymptotics of the functions ip n {z) is 
obtained in the whole complex z plane by solving a matrix Riemann-Hilbert 
problem, following earlier works j5l) |5^] . Very useful is the non- linear steepest 
descent method devised by Deift and Zhou |36| . This set of works Jl2| |5?J 



39 not only provide rigorous and more general solutions for methods and 
ansatzes previously used, but it seem to provide answers also for the models 
previously left unsolved (like the asymptotics of recurrence coefficients R n in 
general cases). The recent and difficult developments will require some time to 
be exploited by physicists. 

Finally, while almost all investigations related to the invariant ensemble of 
random matrices considered a probability distribution of the form of exponen- 



tial of a polynomial like in eq.(l.l), with possible addition of logarithmic terms, 
like the Penner model and the Kontsevich model, there exist probability dis- 
tributions, still invariant under diagonalization by unitary matrices, which lead 
to different forms for the connected correlators Then the classification of 
universality classes perhaps is not yet complete, even in invariant one-matrix 
hermitian ensembles. 

2 Further matrix ensembles and singular value 
decomposition in complex non- hermitian ran- 
dom matrices 

Random matrix models more general than the hermitian one-matrix invariant 
ensemble are often more interesting because of the possibility to describe more 
interesting statistical models. The analytic solution of multi-matrix models 
both in the " perturbative phase" or in different phases is more complex. While 
an ensemble of hermitian random matrices describes triangulations of random 
orientable surfaces, multi-matrix ensembles are suitable to describe models of 
classical statistical mechanics on a random two-dimensional lattice. After the 
breakthrough of the Ising model [|7l[ , Jl3[, it was possible to study random 
walks and loops @ , 0(N) model fffi M g§ @ @ , the Potts model 
fS3[ [105 , surfaces with holes ]73| , [|30|7"a special case of 8-vertex model 
106 , the chiral random matrix which simulate the spontaneously broken 
phase transition of QCD (l0| @ |0§ || |fl). Often it was possible to 



evaluate the critical exponents at phase transition. Multi-matrix models of her- 
mitian matrices are also a good framework for combinatorial problems like the 
four-color theorem |^| [^4| |Q or the enumeration of meanders |H] . Most 
influential, for quantum field theorists were the phase transitions in models of 
random unitary matrices, describing one-plaquette of the lattice formulation of 
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QCD (57) |67|] , the saddle-point solution of the one-matrix ensemble [ fl6[ , 
the Witten conjecture of a master field [104|. 

In this section I shall recall the singular value decomposition, which plays 
a role in the analysis of models with rectangular random matrices and square 
complex non-hermitian matrices. 

Let 4> = (4>ij) *=!. -.« be rectangular random matrix with entries real 

j = i , . . , JVf 

or complex numbers, and joint probability distribution P((f>) invariant under 
<fi — > U 4>V , with U unitary of order N , V unitary of order M 

P(0)# = e -^r 1/(0+0) # / f e -NTrV(4>U) 



V{4>^) = \a 2 {4>H) + \ ai {4>Hf + .. + ^a 2p (^4>Y , a 2p >0 , 
d(f) = Y\ d 4>ij j if is real; or 

i = l,..,N 
3=1, ..,M 

d(j) = Y\ (dRe<ftij dlm<ftij) , if (j) is complex (2-1) 



i=l,..,N 
j = l,..,M 



The hermitian matrices <p'<p and <j)(jy are positive semi-definite, have the same 
non-vanishing eigenvalues fcj = af , where <Ti are the singular values of <f> , 
themselves positive definite. It is straightforward to evaluate the statistics of 
the singular values in the limit N — > oo , M — > oo, while the ratio N/M = L is 
fixed, by the ordinary sa ddl e point analysis. Let us consider first L > 1. The 
probability distribution (2.1) is 

i=l,..,M I i=l,..,M 

where the Jacobian J(iJi) may be evaluated with help from |(3(J 
J(aO = J] (ai) N - M [] if is real, 

i=l,..,M l<i<j<M 

= II (^) 2W " 2M+1 II l^- CT H 2 ' if is complex 

i=l,..,M l<i<j<M 



For the simple case V{<$4>) = \a,2{<ft 4>) + ja^cft (j)) 2 one easily obtains 29 
complex rectangular matrix (f> 

G{z) = lim j-Ttm < ~~7t~7 >= / dt^- , < A < B 

M z — </>'</> J A z — t 



for 
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= ±y/(B-t)(t-A)[^ 



a 4 (A + B) 
8t 



«2 

2t 



(2.2) 



The extrema A, B are given by the usual couple of algebraic equation. Then 
for a 4 = one finds the distribution of singular values for rectangular matrices, 
which is the generalization of Wigner " semi-circle law" 



a 2 y/(B-a*)(a*-A) 



A 



VI- i 



«2 




Returning now to eq.(2.2) for square complex matrices , L = 1, one finds a 
"perturbative" phase for a 2 > — 2-y/aJ with A = where observables correspond 
to resummation of planar graphs, and a " non-perturbative" phase for a 2 < 
— 2-y/aJ where A > , quite similar to the random hermitian case. These results 
were rediscovered by several authors Q 
matrix H and the partition function Z 



MM who introduced the hermitian 



H = 



<j) 
0t 

M 

\\dxi e 



Z = DH e -fr*v{HiX) 



-2/3 V(x- 



I2W-2M+1 



n (^-^) 2 



= 1 l<i<J<M 

The eigenvalues Xi of the " chiral" matrix H are in two to one correspondence 
with the singular values Oi of A : Xi = ±0^. The technique of using the 
auxiliary matrix H in the study of square complex non-hermitian matrix ip was 
developed by Q Q [pll into a powerful method to obtain the distribution 
p(x, y) of complex eigenvalues A^ , see also the similar and simultaneous paper 
|4) . Indeed 



P{x,y) = 



1 N 

J^Yl < *(«-ReAi)5(j/-ImAi) >= 

i=l 

\ d d 1 



■k d z d z* N 
1 d d 1 
tt dz d z* N 



<Tr (JV) log(z - </>)(z* - ft) >-- 
[<Tr {2N) \ogH>-tirN 2 ] 



where now 



H 







4' — z* 

Diagrammatic rules may then be used to evaluate the resolvent G{ri; z, z* 
1 1 n 1 



G(v,z,z*) = — < Tr 



(2iV) 



1 77 



N 



T? - (z*-ft)(z- 



> 
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Finally in terms of the integrated density of eigenvalues of H , Q(p,;z, z*) = 
< Tr(2Ar)0(/^ — H ) > , Feinberg and Zee Q obtain 

As specific examples, Feinberg and Zee consider probability distributions P{<t>, <jy) 



for the complex matrix <p of the form (2.1), invariant under <j) — ► e 



e~ la (jy , where it is natural to expect that the distribution of complex eigenval- 
ues p(x,y) has rotational symmetry p(x, y) — p(r)/(2w), as it indeed happens 
with the Ginibre gaussian ensemble |5(|. For the simple model with only two 
coefficients 0,2 , 04 they find that for 0,2 > — the complex eigenvalues fill 
(non-uniformly) a disk centered at the origin, while for 02 < they fill a 

ring. This is expected from the analogous distribution of the singular values of 
the matrix <p . Next they prove the surprising "sing le r ing theorem" asserting 
that for a generic polynomial potential of the form even in the multiple- 

well cases where the distribution of the singular values of the matrix (j) has 
support on several segments of the positive real axis, the distribution of eigen- 
values of the matrix (f> may only have one ring at most. It seems to me that in 
such cases the assumption of rotational symmetry should be checked. If it turns 
out, by using an explicitly symmetry breaking term to be removed after the 
thermodynamic limit N — > 00, that one obtains different distributions, depen- 
dent on the direction of the symmetry breaking probe, it would be a remarkable 
example of spontaneous symmetry breaking of rotational symmetry. 



3 Conclusion. 

Phase transitions generically occur in the study of ensembles of random ma- 
trices, as the parameters in the joint probability distribution of the random 
variables are varied. They are important for the physics interpretation of sta- 
tistical models and they affect all the main features of random matrix theory 
In invariant one-matrix ensembles recent progress of mathematicians seems to 
solve long standing problems related to multi-cut solutions with no symmetry. 
Phase transitions in ensembles of complex non-hermitian matrices were recently 
explored and it is likely that a richer variety of phase transitions will be discov- 
ered. Random matrix ensembles with a preferential basis, like band matrices 
were studied since the beginning of random matrix theory , see for instance [[sof 
[ jiT) Isij |p9| . Even the simplest cases of tridiagonal matrices with random site 
or random hopping could be analytically solved only for a very limited choice 
of the probability distribution . Well known discrepancies between the moment 
method and numerical methods suggest the presence of phase transitions which 
seem more difficult to analyse than in case of invariant ensembles. 
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